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Abstract 

The usual renormalization procedure for the variational approximation with a trial Gaussian ansatz 
for the A</> 4 model in 3+1 dimensions is re-analysed as a departing framework for the investigation of 
the parameters of the model. The so-called asymmetric phase of the model (where < vac\<j)\vac >^/= 0) 
is considered for the search of privileged values of these parameters (mass and coupling constant) 
and possible conditions they may be expected to satisfy. This also may yield a suitable approach 
for the investigation of the reliability and stability of the approximation. The extremization of the 
renormalized energy density with relation to the renormalized mass, coupling and <j> is done. The 
minimizations of the renormalized energy with relation to the mass and </> provide different expressions 
from the ones obtained by the usual variational procedure for the regularized theory. Sort of “energy 
scale” invariances in expressions for the renormalized mass and coupling constant are found. A different 
view on the restoration of symmetry issue is discussed. The transcendental character of the GAP 
equation may be reduced or even eliminated by placing some variables in the complex plane. 

Key words: Mass, coupling constant, variational method, symmetry restoration, vacuum, non per¬ 
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1 Introduction 


The most developped approach to solve interacting field (and many body) theories is perturbation theory 
which only works well for very small coupling constants as it occurs in Quantum Electrodynamics. Also 
in this approach there is a systematic and direct way of dealing with ultraviolet (UV) divergences, i.e., 
one knows precisely how to renormalize the parameters and to make the theory finite jl], j|, ||. There 
are many motivations for the development of non perturbative methods in Quantum Field and Many 
Body Theories such as for the description of strong interacting systems (with or without spontaneous 
symmetry breaking(s) (SSB)), bound states and phase transitions. One method which has been quite 
extensively investigated is the variational approximation which, with the use of Gaussian wave functional, 
has been showed to be useful in a wide variety of situations. It corresponds to a summation of “cactus” 
type loop diagrams 0, § 0 0. It is equivalent to the Hartree Bogoliubov approach [|j] and also 

to the leading order large N approximation I Q- In this approach the ground state of the system is 
determined by equations for the variational parameters, which are choosen to be a mass and the classical 
expected held characteristic from a SSB state (< p >= <p , which will be referred to as condensate 
[j|]). These equations are derived by the minimization of the (regularized) averaged energy density 
with respect to variational parameters of the trial wavefunctional. The subtraction of the (equations 
of the) theory with < p >= 0 from the ones of the theory in which < p> >=/=■ 0 provides a general 


and consistent elimination of the ultraviolet divergences. Limitations pointed out and discussed in [ 121 
have been rediscussed, leading to extensions and higher order calculations for static and time dependent 
formulations [l3| , 14, 0 m 0 in, 0 1 s 0 0, i 


24, 25 


27, 11]. For the sake of conciseness 


they will not be discussed here although the main ideas of the present work apply to complementary 
approaches. 

It is usually highly desirable to predict the values of the free parameters of a physical theory, such 
as masses and couplings, from the theory itself before comparisons to experimental observations. For 
this there may be possible to predict values, either exact ones or ’’privileged” range. These ’’privileged” 
values may also be associated to the validity of the approximation method or even about applicability of 
the model. The main aim of the present work is to suggest and investigate one reasoning according to 
which values (or range of values) for these parameters could be found. Eventually this may constraint 
them in a more specific theory. Basic ideas are to search for renormalized couplings and masses which 
extremize (minimize/maximize) the renormalized energy density. This procedure can be considered 
as complementary to the renormalization group method [p3], Another procedure will be to consider 
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parameters in the complex plane to introduce auxiliar (imaginary) variables which are to be eliminated. 

The A</ 4 model has been extensively studied for different reasons among which to shed light on non 
perturbative effects in quantum field and many body theories (QFT, QMBT). It corresponds to one of 
the simplest self interacting model whose structure is expected to be (partially) present in several more 
elaborated theories and it presents interesting features [|^, 30, |7], 28, 31, [3|, |33|, 34]. It has also been 
considered for the study of cosmological models [^] and of the Higgs particle in the standard model, for 
example in [^] . It also shares several properties with the linear sigma model (LSM) which is an effective 
model for low energy QCD. Although it strongly seems to possess asymptotic freedom in the asymmetric 
phase the model is “trivial” in the symmetric phase [||, |34|, ^8| . 

In the present work the usual renormalization scheme of the Gaussian approach as carried out, for 
example, in |3(| for the A</ 4 model is used as starting point for further investigation. It is proposed the 
extremization of the renormalized energy density with relation to the renormalized parameters (coupling 
constant and mass). Besides that some variables are placed in the complex plane to search suitable 
(physical) values and eventual conditions for these parameters. The work is organized as follows. In 
the next section the Gaussian approximation is summarized: the GAP equation (transcendental) is 
derived, obtained from the regularized theory (with a cutoff). The renormalization procedure of the 
mass and coupling constant as proposed in |3f| is considered. In sections 3, 4 and 5 values of the 
renormalized mass, condensate and coupling constant which extremize the energy density are searched 
and analysed. They also could yield privileged values of the parameters with which the approximation 
may be more appropriated. In some cases instabilities are found for values of the parameters. In section 6 
a mathematical trick is used to search non transcendetal solutions or/and an expression which constrains 
further the parameters. This is done by allowing some parameters to be complex such that the imaginary 
part must in fact disappear in the end of the calculation. In the last section the results are summarized. 


2 Gaussian approximation for the A0 4 model 

The Lagrangian density for the scalar field </>(x) with bare mass mg and coupling constant A is given by: 

A x ) = \ |^0(x)9 /i </>(x) - m^ 2 (x) - A^ 4 ( X )| (1) 

The theory is quantized in the Schrodinger picture ]37j being the action of the field and momentum 
operators over a state | T [</] > given respectively by: 


<t >|T >= 0|T > 


7r = —iTi-— IT > 


( 2 ) 
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In the static Gaussian approximation at zero temperature the trial ground state wave functional 


is parametrized by the Gaussian: 


V [</>(*)]= Nexp { —J dxdy5(/)(x)G 1 (x, y)^(y)| , (3) 

Where <5</>(x) = 0(x) — </>(x) is the field shifted by the condensate, the point where the wave function is 
centered; the normalization factor is N, the variational parameters are the (classical) expected value of 
the field, <^(x) =< 'P|</i>|'I' >, and the quantum fluctuations represented by the two point function, i.e., 
the width of the Gaussian: G(x,y) =< \D r |0(x)</>(y)|\E r >. In variational calculations the averaged energy 
calculated with \k[0(x)] is to be minimized to obtain the GAP equations. In principle it would yield 
a maximum bound for the ground state (averaged) energy, although ultraviolet divergences make this 
not necessarily reliable. The minimization of the renormalized theory may be useful for this theoretical 
bound of the variational principle. Each of these variational parameters represents one component of the 
scalar field: the expected value in the ground state (’’classical” part) and the two-point Green’s function 
with the mass of the quantum which is decomposed into creation and annihilation operators [§]. 

The average value of the Hamiltonian is calculated and expressed in terms of the variational param¬ 
eters by means of expressions (|2) and (j3j). It is given by: 


n 


= \ \G 4 (x, x) - AG(x, x) + mgG(x, x) + |G 2 (x, x)+ 
+ m o^ 2 ( x ) + (V</>(x)) 2 + f^ 4 (x) + ^<p(x)G(x, x) . 


(4) 


Although in this expression the variational parameters were allowed to have spatial dependence they will 
be assumed to be constant. Variations of the averaged energy density with respect to the variational 
parameters yield the following GAP and condensate equations which define the ground state of the model: 


sn 


0 = -^G 2 (x,y) + + ^(x) 2 (t) 


<5G(x,y) 

-*• 0 = r(x, y)^(y) + ^ 2 (x), 


(**) 


(5) 


Where T(x, y) = — A + [tuq + ^G(x,x)j S(x — y). The Green’s function G may be written from expres¬ 
sions above as: 

G 0 (x,y) =< x| 1 =|y > (6) 

V —A + m z 

where m? is given by the self consistent (transcendental) GAP equation (expression (|j)): 


m 2 = m o + ^ TraceG(x,x,m 2 ) + 


(7) 
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An analogous expression holds for the case in which cj) = 0, i.e., 

/i 2 = m 2 ((f) = 0) = tuq + — TraceG(x , x, /i 2 ). 

Expression (6j) is equivalent to the Feynman Green’s function with time integrated and with sign changed 
in the imaginary part by replacing the self consistent mass by the bare mass ?Bq. The physical masses 
in the different phases may assume different values from each other. The condition of minimum for this 
procedure and its stability was partially investigated in |7j and it corresponds to analysing the second 
order variation of the energy density with respect to the variational parameters. 

From the above expressions it is seen that the non zero solutions for the condensate, cj), can be written 
as: 

4 ? = _ 6 =° - 3G(m 2 ) = 5=1. (8) 

A A 

For G = 0 the tree level value for cj) is obtained in terms of the bare mass. 

The above expression for the Gaussian width @ (and its inverse G 0 1 ) can be calculated in the 
momentum space with a regulator A (cutoff) yielding (for A >> m 2 ): 

a{ "‘ 2 ) = sM a2 - mHn (^)) ’ 

G_1(m2) = 8? (“ J + - x - m ' Ln (Js)) ■ 

where d = 2/ y/e. In the (local) limit of infinite cutoff the average energy and observables diverge and the 
divergences must be eliminated. The renormalization procedure has been performed in three dimensions 
for example in © IH0- 


2.1 Renormalized parameters 


The renormalization procedure of the parameters of the model is done as follows. The energy density of 
the symmetric phase, as well as its GAP equation ([/]), is subtracted from the corresponding expression 
of the asymmetric phase. The GAP equation as defined in expression ([?]) can be rewritten as: 

^ = m 2 + g R (^) 2 + ^ Ln > ( 10 ) 

where the renormalized parameters were defined as: 


2 2 _ 
H =m R 


rriQ + 


AA 2 

167T 2 


1 + 


A 

16 ir 2 


9R = 


Mv 


1 + 


A 

167T 2 


log (f) 


(11) 
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In the first of these expressions m 2 R = g 2 was chosen to produce the usual effective potential f| 36]. It is 
seen from the second of these expressions that in the limit of A —► oo the bare coupling constant would 
go to zero in order to keep gn finite if g is kept constant. This is the “triviality problem. 

The resulting subtracted energy density, TL su b = 'H{4>) — T~C(4> = 0), is re-written in terms of the 
renormalized mass, coupling constant and the mass scale eliminating the cutoff. It is given by: 


H-sub — 


+ 


m 2 TO 1 / 2 2\ 2 1 ( 4 r ( 771,4 

■ (m — g H-— -7T m Ln — T 

4g R \ ** ) 128vr 2 l l g± 


4 i 4 
— m + g 


( 12 ) 


The mass scale g 2 is not a free parameter in fact, it can be considered to be a function of the mass m 2 


and the coupling g R by the GAP expression (10). Other approaches may be of interest for investigating 
the variational method in the Schrodinger picture [3]. In the ground state the parameters (j),m 2 ,g 2 (for 
a given g R ) are related by the GAP and condensate expressions shown above. Any deviation of the 
respective numerical values from the ones related by these expressions induce temporal evolution [||. 

It is possible to verify whether the renormalized GAP equation obtained from the regularized energy 
density given by expression (1C) still is a GAP equation for the energy density given by expression ( 0 ) in 


two ways. The minimization of expression (12) with relation to m 2 is done in the next section. However 


the integration of the GAP equation with relation to m 2 should result in an expression equal to (12) if 
the order of performing renormalization and extracting the ground state does not change results. The 
integral of the GAP equation is given by: 

12 , »» rG"MII dm 2 = 


-g 2 + m 2 +g R 


m 2 / m 

+ s^ Ln V J 


2 2 1 m 1 72 2 1 SR 

= —g m + — + gR(p m + 


m 4 r ( m 2 
-Ln —rr 


m 


(13) 


+ C((j), g 2 


n 2 \g 2 ) 6 g 2 / 

where C(cj), g 2 ) does not depend on m 2 . This expression contains terms very different from the renormal¬ 
ized expression (12). This means either that the minimization of the regularized energy is not equivalent 
to the minimization of the renormalized one or/and that the renormalization procedure has to be im¬ 
proved to make both procedures coincident - if this is possible or desirable. This will be discussed below 
with the minimization of the energy density with respect to (j). 


3 Energy density and renormalized mass 


In this section the renormalized energy density 7i su b is extremized with relation to the renormalized 
(physical) mass: 


dHsub 

dm 


= 0. 


(14) 
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The roots of the resulting expression were calculated by considering that q? is in fact dependent on m? R 
by expression ©■ The resulting expression is given by: 


0 = m 3 

where a* can be given in terms of 


Ln 2 


oi + Ln a 2 + a 3 


(15) 


J = 1- 


9R 

(8tt) : 


= 1 -Gr, 


by: 


ai = — J 2 + 


1 


a 2 = 


Si? 

2 

9R 


327T 2 ’ 
— 1 + J + 


J 2 

(32vr 2 


+ 


128tt 2 


1 + 


2 J 2 
(8vr) 


(16) 


03 = 


32vr 2 


1 + 


9r 

32tt 2 


Expression (ED is not equal to the GAP ([15]) obtained from the minimization of the regularized energy 
density with relation to G(m 2 ). There are therefore five solutions for the renormalized mass m? which 
can be written in the following form: 


where: 


H ± = 


m A = 0, 

m± = n exp^H^ 1 ), 

— o 2 ± a\ — 4aiao 


2ai 


(17) 


(18) 


These solutions for rn^ can be viewed as having corrections for the value of fi due to the self interaction 
through the parameters H± due to the appearance of </> / 0. It is noted that there is a sort of “energy 
scale” invariance in these expressions for rri^ with simultaneous changes in the mass renormalization 
parameter /i. 

The particular case of m 2 = /x 2 , for which the GAP equation is not necessarily valid because cj) —> 0, 
is found for 

a 0 = 0, -> g R = -327T 2 . (19) 

This point may correspond to a restoration of the symmetry. 

The zero mass solutions correspond to a saddle point, they are not minima neither maxima of the 
energy density. If the others solutions are minima is checked via the positiveness of the second derivative: 


d 2 H su b rn 2 


m 


2 Ln ( — ) a\ + <22 ) > 0. 


( 20 ) 


dm? (8vr) 2 

For the derivation of these expressions the complete self consistency of the Gaussian equations was not 
considered. There has been used a truncation on the dependence on /x, i.e., the dependence of Ln(fi/m) 
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on g (self consistency) was considered only for g not very different from m, i.e. g 2 = m 2 + 6 where 
5 « m 2 . Out of this range the above solutions are not expected to be valid. 

In Figures la and lb the solutions of the above equations (m±/g from ©) are shown as a function 
of Gr = <?r/(327t 2 ). All the solutions of figure la, for m + , correspond to stable solutions ( d 2 TL/dm 2 > 0). 
The solutions of figure lb, for m~, are stable for Gr nearly equal or smaller than —1.45 or equal or greater 
than nearly 1.25. The point gR = 0 is not plotted. Values between — 1 < Gr < 0 do not correspond to 
physical stable values of the condensate as it will be shown below, in expression (^J). In the limits of 
gR —► ±oo we obtain analytically that either m = g or m = 0. For the case g —► oo the renormalized 
coupling constant Gr —> 0. While the solution m R in the weak coupling regime can be identified to the 
renormalization point usually considered (for g » m and/or the cutoff going to infinite) there is another 
stable solution m + for which g ~ m + . 


4 The condensate: 0 


The variational equation for the condensate (expression (|j (ii))) is obtained from the regularized energy 
density TC re g- The minimization of the renormalized energy density is done in the following: 

^= 0 . ( 21 ) 

For this derivation the GAP equation provides the dependence of the mass on the condensate, i.e., m 2 ((f>) 
and g 2 = m 2 ((j) = 0) is kept constant. It yields the following expressions: 


6 = 0 , 

6 2 = -^(l + 
9R 


1 

87T 2 


Ln ( — 


m 


( 22 ) 


This last expression may coincide with the expression of 0o obtained from the minimization of the 
regularized energy density (expression (||)) depending on the relation between A and mass scale g as 


it will be shown below. However it is not completely consistent with the GAP equation (1C) which is 
obtained from the minimization of the reguralized energy density with respect to the mass m in the 
asymmetric phase and then renormalized. To make these expressions compatible it would be necessary 
to consider the following alternatives for these expressions: 


g 2 / m 2 R , or g 2 = (g R - 1 )^Ln , (23) 

where m? R is the one of expression ( 0 )- It is not clear whether these identifications are reasonable 
or if they imply a meaningful loss of generality. The limit of gR = 1 does not seem to be reasonable. 
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Therefore the two minimization procedures (of the regularized and the renormalized energy densities with 
respect to the regularized and renormalized parameters respectively) do not seem to yield necessarily the 
same expressions for the parameters in the ground state. Nevertheless it is worth to remember that 
renormalization is performed basically from the regularized GAP equation. 


From the expression (22) the following conditions to obtain real values of cf> can be considered: 


if ■ 9R> 0 ->• Ln 
if : 9R<®—* Ln I 


m 
( m 


V l> 


< -8vr 2 , 


> —8ir z . 


(24) 


The energy density is expected be stable for the condensate values found in expression (22). This 
minimum is verified by calculating the second derivative of the energy density with relation to </>, i.e.: 
d 2 Tt S ub/dp 2 > 0. Its positiveness corresponds to the condition: 


9r[1 + 32 ^) 


> 0. 


(25) 


From this it is seen that for positive coupling constant g r, it can assume any value (from this stability 
criterium) whereas if Qr < 0 one would have to consider gji < —327r 2 . Expressions ( |24| ) and (^) may 
correspond to constraints for the values that the renormalized coupling may assume in order to yield 
stable real ground states. 


Expression (22) can be written as: 


9R(t > 2 = ~m 2 1 + 


1 

8tt2 


Ln 


(26) 


When g = m exp(S-K 2 ) it follows that either f = 0 or gn = 0 in the asymmetric phase of the potential. 
This may correspond to the so called symmetry restoration when the condensate disappears at a partic¬ 
ularly high excitation energy, i.e., the symmetry is restored. A different solution for the particular limit 
of (j) = 0 was found in expression (if) where the energy density is minimum with relation to the mass for 
m 2 = g 2 . 

The above expression for the condensate (^) can be equated to the previous (regularized) one (j|). 
Taking into account the expression of the renormalized coupling constant in terms of the bare one 
(expression (^TJ)) this can be written as: 


A = 


16tt 2 


Ln 


(f) 


-1 + 


2 ( 1 + sh i "(: 


(27) 









If the cutoff is sent to infinite the bare coupling constant assumes different values depending on the ratio 
of p/m. For example, there is a case in which A = 0 if either A —> oo for finite p or: 

— = exp(4ir 2 ), (28) 

being therefore m 2 » / 1 2 . Varying p together with A there may have non zero A solutions. For A/p 
finite, the coupling A may even diverge when: 

— = exp(— 87t 2 ). (29) 

A* 

This is the same point found above (for expression (^) for the possible restoration of the symmetry. 

It is worth emphasizing that it has been assumed, as it usually is, that the minimum of the effective 
potential with relation to the condensate coincides necessarily with its minimum in respect to the physical 
mass m 2 r{ in the regularized theory. 


5 Analysis of the renormalized coupling constant 


Analogously to what was done for the renormalized mass in the preceeding section the extremization of 
the renormalized energy density with respect to the renormalized coupling constant is done in this section. 
Moreover one relevant subject for any approximation method is the understanding of the range of values 
of the parameters of the model (as mass and mainly coupling constants) for which the approximation is 
more appropriated. The extremization is found from: 

&H sub = 
dg R ~ ' 


It is considered, in the following, a truncation of the self consistency of the GAP equations. This is done 
by taking the scale parameter to be close to the mass p 2 = m? + <5, where 5 « m 2 is determined from 
the GAP equation self consistently. From the renormalized GAP equation (expression (|l(^)) it follows 
that: 


5 = 


9R(t> 2 


1 + - 2 & 


16tt 2 


(30) 


Since 6 « m 2 either cj) is large or g r is very large for positive coupling g r. The minimization of the 
renormalized energy yields the following third order algebraic expression: 


(G' r ) 3 + ( G’ r ) 2 (3 + (1 + H)^j + G r (3 + (1 + H)1 + (1 + H)^j + 1 + (1 ^ = 0, (31) 
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where 


H = 


Ln ( 


(87T 2 


s </ _ qr 

Gr ~ (16^2) 


In figures 2a, 2b and 2c the solutions of expression (|3l|) are showed as function of a limited range of 
H, i.e., Ln{m/g). Figures 2b and 2c exhibit the same behavior. It is plotted only the region in which 
the above truncation scheme of the self consistency may be expected to be reliable. The values for gn are 
large and obtained for m 2 ~ g 2 which cannot be simultaneously compatible with the results of Figures 
la and lb (from the minimization of the energy with relation to the mass rn 2 ) - in this limit of <5 << m 2 . 
The point H = —0.001 corresponds to m/g = 0.985 which is not obtained for the larger values of g r 
from figures la and lb. 

These solutions ([H|) may nevertheless correspond to minima or maxima of the solutions of the above 
equation is also verified. For this the positiveness of the second derivative is calculated: 9 ^l ub > 0. All 
the solutions have a negative second derivative corresponding to maxima, instead of minima, of 'H su b{gR)- 


The solutions for the coupling constant of expression (31) depend only on the ratio m/g and not 
on the absolute values of these parameters. This may be seen as a sort of energy scale invariance 
for different physical processes (eventually in different systems) at different energy scales with different 
physical masses. 


5.1 Fixing the energy density 


For the analysis of the system with an energy density given by H su b and a given mass scale, /r, (acceptable) 
values for the renormalized coupling constant and mass may be suggested such that the values remain in 
the physically allowed part of the phase space of the model |27j . This corresponds to fix the renormalized 
mass and energy density (TL su b) and to calculate the resulting physical coupling constants for the process 
involved at a scale g. A third degree algebraic equation is obtained and which can be written as: 


9r 


H 2 


128tt 2 


+ 9r 


H 2 


2 H 2 ' 


x 


+ 9R\- 


H(H + 1) 1 (H 


+ 


- 1 + H‘ 


2 ^-^-- + — = 0 . 


rrr 




(32) 


where x = 1287T 2 and H = Ln(m/g)/(8rr 2 ). This expression also presents a sort of “energy scale” 
invariance for the parameters m/g unless for the term which depends on the total energy density, if Ti su b 
scales differently from rrr. 

In figures 3a, 3b and 3c the solutions of this algebraic equation are shown as functions of H for a 
fixed energy density TL = (100MeF)~ 4 and uir = lOOMeF. The coupling gn can be strong in the region 
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of g ~ m. In particular in the limit of \i = m an unique value is found, it is given by: 

1 


9R = 


/ 'J-^sub 


1 


(33) 


("nF + 1287?. 

For H —► — oo, which is equivalent to g/m —> oo, it follows gu 0 as seen in figure 3a. However solutions 
of figures 3b and 3c do not correspond to gn —> 0, but to a finite (quite strong) value close to 10. 


6 The transcendental character of the GAP equation 

In this section an heuristic trick is used to extract analytical non transcendental solutions from the GAP 
equation or to provide possible further relations among the parameters. Firstly it is considered that the 
the mass scale and renormalized mass may have imaginary parts: 


9 


v 2 = re i9 , 


m 


t = te 


(34) 


where r, s, 0, uj are respectively modulus and phases. With these parametrizations the GAP equation © 
can be written as: 


(rcosO — tcosu) — gR(j) 2 — D(t cosuj Ln(t/r) — (u> — 9)tsinuj )j + 

+i [t sinui — r sinO + D(ujt cosuj — tO cosuj + t sinu Ln(t/r))\ = 0, 


(35) 


where D = . Both the real and the imaginary parts in this expression have basically the same structure 

of the usual GAP equation. It is worth emphasizing that requiring the imaginary part to disappear in 
the GAP equation implies real mass parameters if only one of them is placed in the imaginary plane, i.e., 
either 9 = 0 or uj = 0. In these cases the parametrization in the complex plane is just a trick to reduce 
the transcendental character of the GAP equation (10). 


Requiring the GAP equation (35) to have only real component (this is considered to be a stable 
system) the imaginary part is set to zero. The expression still is quite complicated but the analysis of 
some particular cases will be very useful. For uj = (1 it follows that: 


rsin9 = —D19, 


This may be written as: 


cos9 = \ 1 — 


B 2 9 2 


(36) 


(37) 


In this case the self consistent character of the GAP equation remains strong. The real part of expression 
(35) keeps the same form of expression ( |l0| ) basically with the mass parameters ?n 2 ,/x 2 replaced by r,t. 
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For 6 = 0 (and lo ^ 0) the resulting expressions for the real part of the GAP equation and its 
imaginary part (to be equated to zero) can be obtained from expression (35). They can be written as: 


r — t coslo — gRcf 2 — D(t coslo Ln(t/r ) — Lot sinoo) = 0, 
+i [t sinoo + D(uit coslo + tLn(t/r) sinu ;)] = 0. 


(38) 


It does not provide simpler solutions and therefore they are not shown. The resulting number of free 
parameters is not reduced because although there is one more expression (Qm(GAP)) there also is one 
extra variable (lo). 

Since the phases are auxiliar parameters it is reasonable to assume they are very small without (great) 
loss of generality for the results. The expression for the imaginary part of the GAP equation in the limit 
when sin(9) ~ 9 and sin(u ) ~ lo is given by: 

uotll + D + DLn =6 (r + Dt ). (39) 


This expression can be regarded as fixing the ratio 9/lo. Several particular cases are analyzed below 
although the more interesting case is obtained for lo, 9 non zero and very small. 

( 1 ) Assuming the phases are equal 9 = lo expression (^) reduces to: 

r — t = D t Ln , (40) 

which fixes the ratio r/t or correspondent^ m 2 /g?. This expression is only consistent with the renor¬ 
malized GAP equation for (j) = 0 (which is obtained from the minimization of the regularized energy 
density). Besides that it was mentioned above that, since lo ^ 0 and 9 ^ 0, it is not clear whether g? 
and m 2 remain real although the GAP equation is necessarily real. This may happen because in this 
case the imaginary part of both parameters may cancel with each other to result a real GAP equation 
instead of allowing for independent cancelation. On the other hand each angle (lo or 9) may be set to 
zero separatedly. 

(2) For lo = 0 it follows: 

r ~ — Dt, (41) 


which also fixes the ratio m 2 / g? being a real number only for gR < 0. 

( 3 ) For 9 = 0, expression (39) is computed up to the order of 0(lo 2 ) and it reduces to: 


lo 2 = 2D + 2 + 2 Ln 


(42) 
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where it has been assumed that sinuo ~ lo. In this case it is reasonable to consider cu 2 ~ 0 leading to the 
expression: 


t 


- = exp (—2 — D ). 


(43) 


If lo 7 ^ 0 it will appear in the real part of the GAP equation and therefore the number of free 
parameters in the renormalized equation does not diminish with the new parametrization. Therefore 
lo = 0 would be the only possibly interesting case. This does not happens because of expression ([0]) 
which imposes negative coupling gn < 0. 

The real part of the GAP equation for small angles keeps nearly the form of the original GAP, it can 
be written as: 


t - r + g R (j) 2 + D 


tLn ( - ) — t lo (lo — 9) 


= 0, 


(44) 


where either r or t can be written as a function of the other by means of the constraints of the imaginary 
parts from the expressions (|40|), (|d]) or (^). In this third case the auxiliar parameter lo was not 
eliminated (although 0 = 0). However for very small phases the expression (^4|) reduces to the usual real 
GAP equation (|10|). In this case the real part of the GAP equation is the same as expression (10) written 
as: 


t — r + gn cj) 2 + DtLn (- 

\r 


= 0. 


(45) 


Simultaneously the renormalized energy density must be a real number. The imaginary part due to 
the introduction of parametrization (|34]) has to disappear. However it is easy to notice from expression 
(39) that the resulting expression for the imaginary part of H su b will be quite complicated. Below 
it will be assumed that the phases have small values. This should not impose great limitations in the 
results because they are auxiliar parameters. With this assumption several simplifications occurs because: 
sin(9 ) ~ 9 and sin(uo ) ~ lo. The result for the imaginary part of the energy density, up to first order in 
the phases, will be given by: 


%m(H su b) = uj 
Where 


t(j) 


tr 


+ 2f 2 A_ + 

2 2 g R 647T 2 


+ 




— r 


0, (46) 


A± = 


1 


1 


4 g R 128vr 2 ' 

One of these variables (A + ) can be identified with the solution of fixed 7i su b(LO = 9 = 0) for /r = m given 


by expression (p3|): 


A<=- 


T~^sub 


m 1 


(47) 


fi=m 
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Expression ([If]) still is very complicated and it can also be used to fix the ratio d/u> which can be equated 
to the same ratio obtained from expression (|39|). However this has been written for 7i su b in the form 
written in (ED, which can be written differently by means of the GAP equation for m 2 = m 2 (g 2 ). This 
allows to write an equation of r as a function of t and eliminate one of these variables. The resulting 
identity reads: 

q, gk + sfeMil-r’ 7 ±2 , 4S . 

0 ¥ + * 2 A- + 4 + (1 + D + DLn(i)) 

In this expression the same parameter is used: D = g #/(87r 2 ). This (highly non transcendental) expression 
appears in addition to the usual real part of the GAP equation, expression (fl5|), making a system of two 
algebraic expressions with two variables (r, t). gn is another input/free parameter. 

7 Summary 

A further analysis of the usual renormalization scheme for the variational Gaussian approximation was 
done. The renormalized energy density was extremized with respect to the renormalized mass and 
coupling and to the condensate. Concerning the extremization with respect to the mass, five solutions 
were found, two of which which can correspond to stable vacua in specific ranges of the renormalized 
coupling constant. For this it was considered that the mass scale g is close to the physical mass. A sort 
of “energy scale” invariant algebraic expression was found in this calculation. In other words, changes 
in the renormalized (physical) mass m 2 with corresponding change in the renormalization mass scale 
parameter g 2 yield the same solutions. 

Values for i j>, in the vacuum, were also found by minimizing the renormalized energy density with 
relation to it. The resulting expression is not completely consistent with the renormalized GAP equation 
unless the expression (Q) is modified such that g 2 / m? R —► 0. From this expression it was pointed out 
that either the “condensate” or gn disappears when the mass scale (introduced in the renormalization 
procedure) assumes the value 

(ft (^g = m exp(8ir 2 )^ = 0 . 

This can be seen as a restoration of the spontaneous symmetry breaking. With this value for g, the bare 
coupling A may also diverge for A/g finite, as shown in expression (|27|). 

Particular values of the renormalized coupling constant which extremize the energy density were also 
found. The coupling constant may constraint the values of the renormalized mass which yield maxima of 
the energy density (effective potential). The extremization of the effective potential with respect to the 
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coupling constant was also performed in the limiting case that the mass scale g is close to the physical 
mass. This is a way of truncating the self-consistency of the approximation. Another kind of “energy 
scale” invariant expression was also obtained. No minima of the energy density with relation to gR 
was found (considered without the whole renormalization group equations) within the truncation scheme 
which was adopted. The renormalized energy was also fixed to provide specific values for the coupling 
constant as a function of the energy density and the mass, which disappears in the limit of g = m. 

The masses were allowed to assume complex values to search non-transcendental solutions for the 
GAP equation and other relations among the parameters reducing the number of free parameters. The 
imaginary part is required to be zero at the end of the calculation producing another expression which 
relates the mass, coupling and the renormalization scale parameter. This parametrization for the imag¬ 
inary part may lead to new relation between the parameters reducing the number of free variables. 
These imaginary parameters may be required to be very small (sin(uj) ~ u> or sin(8 ) ~ 8). The same 
parametrization is applied to the energy density which also must be a real number. The number of free 
parameters ( vn 2 or g 2 , and gR ) is reduced and non transcendental solutions may result such as that of 
expression (|d]). 

The ground state in the framework of the variational approximation is found by the minimization 
with respect to the two point function G(x, y,m 2 ) (which is a function of the physical mass m 2 or 
/ l 2 = m 2 ((j) = 0)) and to the condensate 0 - they are the variational parameters (given in expressions 
(H)). Although they are regarded initially as independent variables, the GAP equation (for ground state) 
relates them. While the GAP equation is used for the renormalization of the bare parameters in the 
vacuum, expression ( |22| ) was calculated from renormalized expressions for the ground state. However 
there is nothing really defined about the behavior of the renormalized parameters in excited states. It 
was shown with sections 2.1, 3 and 4 that the minimizations of the renormalized energy with relation 
to the mass and cj) yield different ground state (GAP) expressions from the ones obtained by the usual 
variational procedure for the regularized theory. This may have several meanings. It may not be evident 
whether these variational parameters are really or completely suitable as independent parameters for the 
Gaussian approximation and extensions (or leading order large N, Hartree Bogoliubov) or even in the 
exact ground state , i.e., the energy must be minimum with respect to particular combination(s) of these 
(or other) (physical?) variables. Notwithstanding the minimization of the regularized energy may not 
be equivalent to the minimization of the renormalized one because in the regularized theory there still 
are other (bare) parameters which are eliminated in the renormalization procedure corresponding to a 
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sort of ’’hidden dependences” among them. It may also be that the renormalization procedure has to 
be improved such as to make both ways of obtaining the ground state expressions equivalent. In this 
case the renormalization procedure would be allowed to be done at any moment independently of the 
order of the the variation, renormalization and extraction of observables within a certain nonperturbative 
approach. 
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Figure captions 


Figure la - First solution of expression (|l^) - a new GAP equation - for the ratio of the renormalized 
mass to the mass scale /r as a function of g R /(8ir 2 ). 

Figure lb - Second solution of expression 0) - a new GAP equation - for the ratio of the renormalized 
mass to the mass scale n as a function of g R /(8ir 2 ). 


Figure 2a - First solution of expression (31) for the renormalized coupling constant - from the min¬ 
imization of the energy density with relation to the renormalized coupling constant - as a function of 
H = Ln(m R /n)/(8TT 2 ). 


Figure 2b - Second solution of expression ([3l]) for the renormalized coupling constant - from the 
minimization of the energy density with relation to the renormalized coupling constant - as a function of 
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H = Ln(m R /fi) / (8n 2 ). 


Figure 2c - Third solution of expression (31) for the renormalized coupling constant - from the min¬ 
imization of the energy density with relation to the renormalized coupling constant - as a function of 
H = Ln(m R /n)/(8Tr 2 ). 


Figure 3a - First solution for the renormalized coupling constant of expression (32) - fixing 7i = 
(lOOMeF) 4 and m R = lOOMeV - as a function of H = Ln(m R /fi)/(8ir 2 ). 


Figure 3b - Second solution for the renormalized coupling constant of expression (32) - fixing Tt = 
(lOOMeF) 4 and m R = 100 MeV - as a function of H = Ln{m R /fx)/{ 87t 2 ). 


Figure 3c - Third solution for the renormalized coupling constant of expression (|32|) - fixing TL = 
(lOOMeF) 4 and m R = lOOMeF - as a function of H = Ln(m R /n)/(8 tt 2 ). 
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